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Abstract. - We propose that the degree of sweehng of the 30nm chromatin fiber (a "mea- 
sure" of its transcriptional activity) is mainly determined by the short-ranged electrostatical 
interaction between different sections of the "folded" DNA chain. These sections constitute 
only a small fraction of the chain and they are located close to the entry-exit points of the 
DNA chain at the nucleosome core particles. We present a model that allows to estimate the 
degree of swelling of chromatin fibers as a function of salt concentration, charge density of the 
strands etc. Different mechanisms by which the state of chromatin can be controlled in vitro 
and in vivo are discussed. 



DNA in eucaryotic cells is organized within a protein-DNA complex known as chromatin. 
In this way plant and animal genomes are compactified into volumes whose linear dimensions 
are many orders of magnitude smaller than their contour lengths. For instance, the human 
genom is made up of billions of base pairs (bp) corresponding to about one meter of DNA 
chains. These highly-charged and hard-to-bend polymers are condensed into complexes that 
have a characteristic size of a micron and therefore fit into the cell nucleus. At the same 
time it is of vital importance that a fraction of the genetic code stored within these tight 
complexes is accessible: For instance, gene regulatory proteins bind to specific sequences, the 
rather bulky protein complex RNA polymerase needs to gain access to a whole gene during 
its transcription etc. How DNA is "folded" within chromatin and how it can be "unfolded" 
for, say, transcription purposes is still poorly understood. 

The primary structure of chromatin is known in great detail from x-ray studies |]. The 
basic unit is the nucleosome consisting of the core particle and the linker DNA (of typical 
length 60bp) that connects to the neighboring core particle. The resulting structure is a 
beads-on-chain necklace (lOnm fiber). The core particle consists of 147bp DNA wrapped in 
1- and 3/4 left-helical turns around a globular octamer of cationic proteins (core histones) - 
forming a squat cylinder with a radius of about bnm and height of about Qnm. The higher- 
order secondary and tertiary structures on scales from lOnm up to a micron are still a matter 
of controversy It is well-known that the bead-on-a-chain folds into a thicker fiber with 
a diameter of roughly 30n?7i but it is still not clear how the neighboring beads are arranged 
in this fiber with respect to each other and where the linker is located. In the solenoid 
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models [| it is assumed that the chain of nucleosomes forms a helical structure with the 
linker DNA being bent whereas the zig-zag- or crossed-linker models Q posits straight linkers 
that connect core particles that are located at opposite sites of the fiber. An example of that 
kind of structure, the two-dimensional zig-zag fiber, is shown in Fig. 1(a). In general, the fiber 
is three-dimensional. 

The difficulties in determining the structure of the 30nm fibers are due to the lack of 
reliable experimental methods. Electron cryomicroscopy allows to visulize fibers in vitro but 
the fiber geometry can only be detected at ionic strengths well beyond physiological conditions 
where the structures are much more open jsj . Fibers at such low ionic strengths show straight 
linkers, a fact that supports the second class of the above mentioned chromatin models - at 
least for fibers at these unphysiological conditions. At higher ionic strengths, however, the 
fiber is so dense that their internal structure remains obscure. 

A different approach allows to obtain structural information via an indirect method: the 
stretching of chromatin fibers with the help of micromanipulation devices as it was achieved 
recently. Cui and Bustamante found that chromatin fibers are extremely soft with respect 
to their extension (compared to naked DNA). In principle, it should be possible to compare 
the obtained experimental data (namely the force-extension curves) with the existing models 
and to determine which is the appropriate model. In the meantime the crossed-linker model 
was indeed reevaluated in view of the stretching experiment. Katritch et al. |^] performed 
computer simulations showing that this model is capable of reproducing experimental data 
quite well for several sets of values of the "free" parameters. Schiessel et al. Q gave an 
analytical treatment and introduced optimization criteria (see below) that led to reasonable 
predictions of the mechanical properties of the chromatin fiber without the use of any adjustable 
parameters. These studies do not invalidate the solenoid model but demonstrate at least that 
the crossed-linker model is successful in describing some characteristic features of the 30nm 
fiber. 

The main idea of the crossed-linker model as it was introduced by Woodcock and cowork- 
ers [Q is that the three-dimensional structure of the 30nm chromatin fiber is determined 
by two quantities only: the entry-exit angle a of the linker DNA at the core particle (c/. 
Fig. 1(a)) and the linker length b. The entry-exit angle a is determined by local properties of 
the nucleosome core particle and it seems reasonable to assume that this quantity is constant 
throughout the fiber as long as there are homogeneous conditions. The second quantity, the 
linker length 6, determines the rotational setting of neighboring nucleosomes. This follows 
from the fact that there are specific binding sites on the histone octamer where the DNA 
adsorbes with its minor groove facing inwards. Taking into account that the helical repeat 
length is « lObp it is clear that increasing b by one bp leads to rotation of a nucleosome by 
K, 36° with respect to its neighbor. As a consequence, changing b by multiples of w lObp does 
not affect the rotational setting, i.e., the dihedral angle is periodic function of the linker 
length. 

It was found that chromatin fibers have the tendency to show a preferential quantization 
with respect of their linker length and that this quantization is set by the helical twist |9| . This 
indicates that the dihedral angle (j) shows a preferential value and thus the relative rotational 
setting is more or less constant throughout the fiber. Therefore it seems to be a reasonable 
approximation to consider both, a and (j) to be constant throughout the fiber (together with 
b). The geometrical properties of the fiber, i.e., its radius R and the linear density A of 
nucleosomes along the fiber, are then a function of the two angles a and (j) only (and of & in a 
trivial sense). In Ref. j|] we termed this geometrical model therefore the two- angle model. A 
specific example of such a two-angle fiber with = 180 is the zig-zag fiber that we depicted 
in Fig. 1(a). 
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Fig. 1 - (a) Schematic view of a section of the 30nm chromatin fiber. For simphcity, the fiber is 
shown as a two-dimensional zig-zag fiber. The histone octamers are displayed as black circular disks. 
The DNA (white) consists of wrapped parts, "stem" sections and sections that link neighboring 
nucleosomes. (b) Enlarged view of the stem region. As shown schematically the stem section of the 
DNA interacts with cationic N-tails of the core histones and with the linker histone HI (see text for 
details). 



In our study of the geometrical and mechanical properties of the two-angle model we 
suggested that the values of a and 4> of the " native" 30nm chromatin fiber in its " silenced" 
form (no transcriptional activity) are chosen in such a way that the three-dimensional density 
of the fiber (defined by K/ttE?) and its ^'accessibility" are optimized. Accessibility means 
here how rapid the "silenced" fiber opens up when its geometrical properties (characterized 
by the two angles (j) and a) are changed for the purpose of, say, transcription. As long as 
the positions of the nucleosomes at the DNA are fixed (as assumed here) <j> and b are always 
constant and hence it is the angle a with which the fiber geometry, i.e., its degree of swelling, 
is controlled. High accessibility means that the fiber opens up strongly with an increasing 
value of a, leading to a strong reduction of the nucleosome line density. Hence our measure 
for accessibility, as introduced in Ref. is the quantity —dA/da. 

The purpose of this paper is to study how the entry-exit angle of the DNA at the nucle- 
oseomes can be controlled via electrostatics. It can been seen in the cryo EM studies |^ that 
the fibers open up and become therefore more accessible when the ionic strength is reduced 
and that this opening is directly linked to an increase in a. It was suggested that via other 
mechanisms (for instance, the acetylation of so-called histone tails ||ic| ] , as explained in more 
detail below) the angle a and therefore the degree of swelling can be changed for a given sec- 
tion of the fiber and that this constitutes a biochemical means to control the transcriptional 
activity of genes. 

Whereas the x-ray studies of the core particle [Q allow a detailed knowledge of the wrapped 
part of DNA it does not give insight into the conformational properties of the entering and 
exiting strands since the core particles were constituted from 146bp DNA only. One has 
therefore to refer to the electron cryomicrographs. In these micrographs it can be seen clearly 
that lOnm stretches of the entering and exiting DNA strands are glued together forming a 
unique "stem motif" [|j (c/. also Fig. 1(a)). The glueing of the two equally charged chains 



4 



EUROPHYSICS LETTERS 




Fig. 2 - Idealized model for the entry-exit region of the DNA at the nucleosome. h [x) and —h {x) 
describe the conformations of the entering and exiting DNA sections, respectively. The asymptotic 
slope at large x determines the entry-exit angle a and, in turn, the overall conformation of the 
chromatin fiber (as depicted in Fig. 1(a)). 

is accomplished - amongst other things - via a special cationic protein, the so-called linker 
histone HI as shown schematically in Fig. 1(b). 

At physiological concentrations (roughly lOOmM salt corresponding to a screening length 
of about lOA) the electrostatics is essentially short-ranged (note that the diameter of DNA 
is 20A). It seems therefore reasonable to assume that a is set within the small region where 
the two linker DNA are in close contact, i.e., in the stem region |jll[. This value of a in 
turn controls the large-scale secondary structure of chromatin, the 30nm-fiber. To mimic 
this situation we assume in the following two parallel DNA strands that are hold together 
tightly at y = for a; < and that are free for a: > 0, cf. Fig. 2. Because of their mutual 
electrostatic repulsion the two strands bend away from each other. When the two strands are 
far enough from each other their interaction is screened so that they asymptotically approach 
straight lines (neglecting thermal fluctuations as it is appropriate for the length scales under 
consideration). The two asymptotes define the opening angle a as indicated in Fig. 2. It is 
the purpose of this study to determine this angle and its dependence on different physical 
parameters [ p^ . 

We describe the conformation of the upper DNA chain by the height function h [x) . By 
symmetry the position of the lower strand is then given by —h[x). From the above given 
considerations follow two boundary conditions at a; = 0: 

h (0) = h' (0) = (1) 

The entry-exit angle a is related to the slope of h (x) at infinity as follows: 

tan(a/2) = /i'(oo) (2) 

We model the two DNA chains as semifiexible polymers with persistence length Ip and line- 
charge density A immersed in a salt solution characterized by the Bjerrum length Ib = /eksT 
and the Debye screening length = {%'kCsIb)~^^'^ (e dielectric constant of the solvent, T 
temperature, Ub Boltzmann constant and salt concentration). We assume that the charges 
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on the DNA strands interact via a screened electrostatic potential. The free energy of the 
system is then given by 



F 



dx 







dx"^ 



2lB\^Ka{2Kh {x)) 



(3) 



P—-2lB\^KKi{2Kh) = Q (4) 



The first term in the integral accounts for the bending of the two DNA strands and the second 
one describes the interaction between the two chains [Kq {x) being the Qth order modified 
Bessel function). Here the interaction of a given charge on one chain with all the charges on 
the other chain is approximated by the interaction of this charge with a straight chain at the 
distance 2h. The integral, eq. ^, is a good approximation as long as the value of a that follows 
from its minimization is sufficiently small, i.e., as long as h' (x) ^ 1 for all x. We furthermore 
neglect here the interaction between charges on the same chain. The conformation of the 
upper chain, h {x) , is then the solution of the corresponding Euler-Lagrange equation (using 
Kl,{x)^-K^ (x)) 

, d-^h 
"d^ 

together with four boundary conditions: Two are given at the origin, eq. |l|, and two follow 
from the condition of straight "linkers" at infinity: 

h" (oo) = h'" (oo) = (5) 

Defining h = 2Kh and introducing the dimcnsionlcss quantity x ~ ^ x, eq. ^ 

can be rewritten as 

^-K,(h)=0 (6) 

It follows immediately 



We denote the dimensionless asymptotic slope by co — dh/dx 
that tan (a/2) is given by 



.an(a/2).V(co,.-i(^) (i) (7) 

Equation ^ can be solved asymptotically for small value of h where Ki (jij ~ \/h. In 

leading order the asymptotic solution is given by h [x) ~ /iq {x) = x^ ^ (— Ini). This solution 
fulfills the boundary conditions at the origin, eq. The boundary conditions at infinity, eq. ^ 
are only marginally affected by the solution for small values of x as can be seen by adding terms 
like ax^ and h5? . We can use this approximate solution to give a rough estimate for cq. (x) 
reaches its maximal slope 0. 56 at a; ~ 0. 41. For larger values of x the approximate solution 
(x) becomes more and more unreliable. This asymptotic analysis indicates that cq is of the 
order one fl^ . The asymptotic solution allows to estimate the free energy of the configuration 
since most of the bending and electrostatic contributions are localized close to the origin. 
Replacing h[x) in eq. ^ by h^ix) = {"fx)^ \/ — In {'^x)/2k with 7 = (^41b\^ /Ip)^^^ and 
integrating from x = Oto x ^ 7"^ (neglecting logarithmic contributions) we find that both, the 
bending contribution and the electrostatic contribution scale as F/ksT « Ig^X^^'^lp'^n"^^'^. 

Our calculation does not account for miramolecular electrostatic contributions. These 
can be included by adding the electrostatic persistence length Iqsf — Ib^^ l^i^'^ (the Odijk- 
Skolnick-Fixman length, cf. Ref. ||l^) to the bare persistence length, i.e., by replacing Ip 
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by Ip + losF- This leads to tan(Q;/2) c>c {Iosf/ ('p + Iqsf))^^^- Note that for smaU salt 
concentrations (large k~^) Iqsf becomes increasingly important so that the angle approaches 
a universal value independent of k^^, A and Is - at least as long as the above stated ap- 
proximations are still applicable in this limit. For the salt concentrations under consideration 
losF <C so that we disregard intramolecular contributions in the following. 

We apply now our results to the problem of how the geometry of the chromatin fiber is 
controlled in vitro and, on a more tentative level, in vivo. The in wifro-experiments show 
that chromatin fibers "open up" with a decreasing salt concentration. From the electron 
cryomicrographs it was found that a « 85° for — 5mAI and a ~ 45° for Cg — 15mM [5|. 
Furthermore, from electron cryotomography it was estimated that a « 35° for Cg = 80mM |5|. 

We expect from eq. |^ that a ~ 2arctan (^Ccs ^^^^ with C being a constant. Let us take the 
angle at the highest salt concentration, = SOmM, as the reference value. From this follows 
C = 0.94. With this value of C we predict a « 51° for Cs = 15mM and a « 64° for 
Cs = 80mA/. Whereas the predicted value of a at the intermediate ionic strength is close to 
the experimental one, the value for low salt concentrations is noticably too low. However, we 
expect the theory to break down for such large values of a. 

How can the degree of swelling of the chromatin fiber be controlled in vivo! Under the 
assumption that the above mentioned geometry is valid the only parameter that might be 
under biochemical control is the linear charge density A. It is known that there are flexible 
cationic tails, namely the C-tail of the linker histone as well as the lysine-rich N-tails of the 
core histones involved in determining the degree of swelling. In Fig. 1(b) we give a tentative 
picture of the conformation of two N-tails that protrude from the histone core (we assume 
here that these are the tails of the two so-called H3 core histones that are located close to 
the entry-exit point). It is known that if either of these components is missing the fiber does 
not fold appropriately (cf. Ref. |0 and references therein). As indicated in the Figure the 
tails might form a complex with the entering and exiting linker DNA in such a way that they 
effectively reduce its linear charge density A. It is known that transcriptionally active regions 
in chromatin show an acetylation of the core histone tails (i.e., the cationic groups of the 
lysines are neutralized) . In our tentative picture this acetylation mechanism would increase A 
and according to eq. ^this would lead to an opening of the entry-exit angle a. The acetylation 
might therefore be the first step in the decondensation of a stretch of the chromatin fiber that 
needs to be accessed for transcription Jist . 

Further steps might then involve the loss of the linker histones leading to a ~ tt, i.e., a bead- 
on- string filament. After the loss of the linker histones the nucleosonies might even become 
mobile as it was reported in Ref. ||l^. We recently suggested that this "nucleosome sliding" 
results from DNA "reptating" around the histone core with the help of intranucleosomal 
loops Loop formation might also be crucial for the actual transcription itself: It was 

suggested that RNA polymerase might elongate through the nucleosome by passing it in a 
loop |pT|] . However, up to now a detailed knowledge of how transcription through chromatin 
is possible is still unclear and a matter of current research. 

* * * 



I wish to thank R. Bruinsma, A. Johner and G. Migliorini for useful discussions. 



REFERENCES 



Helmut Schiessel: How short-ranged electrostatics etc. 



7 



[I] LuGER K., Mader a. W., Richmond R. K., Sargent D. F. and T. J. Richmond, Nature, 
389 (1997) 251. 

[2] Woodcock C. L. and Dimitrov S., Curr. Opim. Genet. Dev., 11 (2001) 130. 
[3] Finch J. T. and Klug A., Proc. Natl. Acad. Sa. USA, 73 (1976) 1897. 

[4] Woodcock C. L., Grigoryev S. A., Horowitz R. A., and Whitaker N., Proc. Natl. Acad. 

Set. USA, 90 (1993) 9021. 
[5] Bednar J., Horowitz R. A., Grigoryev S. A., Carruthers L. M., Hansen J. C., Koster 

A. J. and Woodcock C. L., Proc. Natl. Acad. Sci. USA, 95 (1998) 14173. 
[6] Cui Y. and Bustamante C., Proc. Natl. Acad. Set. USA, 97 (2000) 127. 
[7] Katritch v., Bustamante C. and Olson W. K., J. Mol. Biol, 295 (2000) 29. 
[8] Schiessel H., Gelbart W. M. and Bruinsma R., Biophys. J., 80 (2001) 1940. 
[9] Widom J., Proc. Natl. Acad. Sci. USA, 89 (1992) 1095. 

[10] van Holde K. and Zlatanova J., Proc. Natl. Acad. Sci. USA, 93 (1996) 10548. 

[II] We assume here that internucleosomal interaction can be neglected. There is evidence from the 
stretching experiments [^] that the nucleosomes are "in contact" at a physiological ionic strength 
and that there is a short-range attraction between them. With the swelhng of the fiber at lower 
ionic strength the internucleosomal interaction becomes less important and might be negligible 
(cf. our discussion of the internucleosomal interaction in Ref. |^). 

[12] A similar problem is that of a charged chain that is deflected by a rod perpendicular to it. This 
mimics a situation occuring in semidilute polyelectrolyte solutions and is discussed in Ref. ]l^ , 
p 39. 

[13] Barrat J.-L. and Joanny J.-F., Adv. Chem. Phys., 94 (1996) 1. 

[14] Note that our analysis is based on the assumption of a vanishing distance between the two 
strands for a; < 0. This leads to a logarithmic singularity of the curvature at a; = 0, i.e., to an 
infinite torque at the origin. The structure could only be stabilized for an infinite adsorption 
energy per length between the two strands for a; < 0. A more realistic case is h (x) = e <Si 
for a; < 0. In this case the following approximate solution can be constructed: a;*/24£ — 
/ (3a;o\/— In io) -f x^ \J {— In So) + e for x < xq and Jiq (x) -|- e for xo < x <^ 1; here aio is the 
solution of Xq^J —\n io = £■ These two functions cross over smoothly at a; ~ io and obey Eq. |^ for 
X <^ X(i and for a; ^ io, respectively. This indicates that the solution is only marginally affected 
by the value of e as long as e <C and that our original assumption, Eq. 0, is reasonable. 

[15] The processes that are involved in the acetylation and deacetylation might be quite speciflc and 
involved as, for instance, discussed in Ref. Um. The histone tail modifications might serve specific 
functions via the modiflcation of their secondary structure that in turn modifies their interaction 
with certain proteins jl7|. Recently there are even attempts to decipher a specific "language" of 
covalent histone modifications |l^]. It might be that such specific processes act in concert with 
the more basic charge neutralization principle discussed here. 

[16] Grunstein M., Nature, 389 (1997) 349. 

[17] Hansen J. C., Tse C. and Wolffe A. P., Biochemistry, 37 (1998) 17637. 
[18] Strahl B. D. and Allis D., Nature, 403 (2000) 41. 

[19] Pennings S., Meersseman G. and Bradbury E. M., Proc. Natl. Acad. Sci. USA, 91 (1994) 
10275. 

[20] Schiessel H., Widom J., Bruinsma R. F. and Gelbart W. M., Phys. Rev. Lett, 86 (2001) 
4414. 

[21] Studitsky V. M., Clark D. J. and Felsenfeld G., Cell, 71 (1992) 371. 



